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A thirty-four day record of near-bottom temperature and horizontal 
velocity spanning the lower third of the water column from within 
Monterey Canyon was examined. The observed internal tide is highly non- 
linear with kinetic energy dispersed among numerous overtides near the 
bottom and more concentrated in the primary semi-diurnal constituent 
(M2) higher in the water column. The bottom currents and temperature 
vary in strength over the record period, taking on the characteristics 
of an internal bore at times with large up-canyon accelerations 
accompanied by rapid temperature drops. The bores are nearly phase 
locked to the sea level variations and arrive at the measurement site 
~8.6 h after high tide in Monterey, CA. They are evident in the 
velocity records up to at least 35-m above the bottom and may be 
accompanied by high frequency pulses that extend higher. The variation 
is not caused by direct forcing from the barotropic tidal range as the 
strongest bores do not exclusively occur during either the spring or 
neap phase of the barotropic tide. Speculation on the cause for the 
temporal variation centers on changes in mid-water stratification 
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I. INTRODUCTION 


roe MOTIVATION 

a le Background 

Internal waves are prevalent throughout the world’s 
oceans. An internal wave is a gravity wave within a fluid 
of variable density. Almost any measurement taken in the 
ocean will contain signals of internal waves (Munk 1981). 
Internal waves can even have ae surface signature of 
alternating bands of smooth and rough seas given the correct 
conditions (Munk 1981). Numerous studies have been 
conducted to try and link internal waves to other ocean 
processes. 

Internal waves are thought to play a major role in the 
distribution of energy in the ocean (Munk 1997). An 
empirical energy spectrum (GM79) of internal waves was 
proposed in 1979 (Garrett and Munk 1979) for the open ocean. 
This spectrum has been used as a benchmark for comparison 
with measurements taken in the various parts of the ocean. 
Internal wave energy has been found to exceed the GM79 
spectrum in coastal regions. Huthnance (1989) speculated 
that the GM79 spectrum 1s valid with one localized 
deviations near and on the continental shelf. Eriksen 


(1998) has recorded increased energy levels compared to GM79 


at frequencies that correspond to the local critical 
frequency for sloping terrain. 

Internal waves in submarine canyons exert effects on 
the sediment transport, biology, and even on the water 
properties of the continental shelf farther up-canyon. 
Submarine canyons act like funnels moving water onto the 
continental shelf (Huthnance 1989, Lafuente 1999). Most 
internal waves will be "trapped" and funneled farther up- 
canyon by the walls that are Supercritical for all but the 
shortest period internal waves, while the floor is generally 
subcritical (Summarized in Petruncio et al. 1998). When the 
canyon’s bottom slope happens to be near the critical slope 
for one one the major tidal constituents, bottom 
intensification of the tidal currents, or even an internal 
tidal bore, may be found (Thorpe 1992, Ostrovsky and 
Stepanyants 1989). Petruncio et al. (1998) found the 
Monterey Canyon’s bottom slope to be near critical for the 
M2 tidal frequency, causing intensified bottom tidal 
currents. Shea and Broenkow (1982) observed large magnitude 
internal tides in the canyon and what may have been an 
internal tidal bore near the canyon head. Cacchione (1999) 
has pointed out that the slope of many submarine canyons 1s 
close to being critical for the M2 tide. As such, the 
Monterey Canyon may be representative in its effects on the 


internal wave field. 


2. Experiment Setup 

Submarine canyons are an under-sampled region compared 
to more 2-dimensional areas of the continental shelf. The 
steep topography makes it very difficult to obtain 
oceanographic measurements via moorings or cable-lowered 
instruments. Little is known about internal wave spatial 
scales in canyons where interaction with the steep 
topography is likely to facilitate more rapid transfer of 
energy toward smaller scales and turbulence production. 
Because of their lateral and vertical heterogeneity, canyons 
are difficult to characterize with a few point observations. 
One way to overcome this problem is to conduct profile 
surveys that would provide a broad spatial coverage. A 
National Science Foundation (NSF) funded experiment was 
conducted in the Monterey Canyon by Leslie Rosenfeld (Naval 
Postgraduate School) and Eric Kunze (Univ. of Washington) 
with the following objectives taken from the grant proposal 
for the experiment. 

e Determine the finescale shape and level of the 
vertical wavenumber spectra rex vertical 
displacement and horizontal velocity as compared 
with the open ocean GM79 model. 

e Determine the spatial scales and heterogeneity of 
the internal wave field as a function of vertical, 


along and cross-canyon position. A specific issue 


of interest is how much energy negotiates the sharp 


bends in the canyon versus being back reflected. 


e Determine the anisotropy of the internal wave field 
to evaluate where the dominant semidiurnal tide is 


best described as propagating and where as standing. 


e Intercompare dissipation rates inferred from the up- 
canyon energy-flux divergence, and the finescale 
shear and strain spectra using the Gregg-Henyey- 
Polzin sealing. 

The experiment collected numerous measurements taken 
both by moored and shipboard instruments. For canyon depths 
between 500 and 1600 m, full-water-depth velocity, density 
and dissipation profiles were made using Sippican velocity 
and CTD oprofilers and lLueck’s expendable dissipation 
profilers (which failed to meatracl useful data). In canyon 
depths shallower than 300 m, over the canyon rim and on the 
adjacent shelf region, fine- and microstructure measurements 
were collected by Michael Gregg (Univ. of Washington) in a 
closely associated experiment. In an See to bridge the 
deep and shallow experimental domains, moorings with upward 
facing acoustic Doppler current profilers (ADCP’s) were 
placed between 300-m and 400-m depth to give a month-long 
time series of the lower water column currents. Two 12-h 
time series of CTD casts were made over the canyon near 


400-m depth, and two near 800-m depth. This thesis will 


concentrate on the data collected by CTD’s and moored 


instruments. 


B. INTERNAL WAVE THEORY 

A Linear Internal Wave Theory 

Linear wave theory is a very well understood and widely 
applicable wave theory. Internal gravity waves. are 
described by linear theory as easily as surface gravity 
waves. Internal waves can be formed at the boundary between 
two fluids of different densities or in a continuously 


stratified fluid. 


a) Dispersion Relation 

The dispersion relation defines how a wave type 
behaves. More specifically, the dispersion relation relates 
the angular frequency and the wave number of any given wave. 
This relation is different for every type of wave and the 
dispersion relation for internal gravity waves reveals 
unique characteristics about them. 

Starting with the basic conservation equations 
from Kundu (1990) and including the Boussinesq approximation 


with the inviscid fluid approximation yields: 


Ze, aa 2 (1) 
ot Py Ox 
dv 1 dp 
ee — ae (2 
ot es Py Oy : 


dara apace oa (ay) 
ot oz es 

ox oy az 

PP 6 (5) 
Dt 


Equations (1) through (5) contain five unknowns: 
horizontal velocity (u,v), vertical velocity (w), density 


(0), and pressure (p). Equation (5) 1s a result of the 


assumption that water 1S an incompressible fluid. 
Assuming that the system is starting at a steady 
state of no motion, a hydrostatic balance will exist and the 


vertical acceleration will be initially zero. 


=o=-—2_& (6) 


In order to for waves to appear, motion must 
develop that will disturb the initial pressure and density 
fields. We can describe the motion effects on density and 
pressure in terms of a mean quantity and a fluctuating 
quantity. 

p= p(z)+ p(x, y,z,t) (7) 

p = p(z)+ p(x, y,2,0) (8) 


Now we must use equations (7) and (8) to describe 


a water particle's motion due to the now present wave 


action. This is done by substituting equations (7) and (8) 
ita © the equations of motion (1-5). Because we 
assumed that there wasS some previously existing steady 
State that is now being perturbed, the gradients of the mean 
quantities will be dropped since they only describe the 
state of no motion. Since we are dealing with linear wave 


theory, all of the non-linear terms are also neglected. 





a) tA 
ou _¢,___| (9) 
ot Dyeo 
dv 1 op’ 
—+ fu =-—-— (10) 
ot Py oy 
i i Gey 
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frequency or buoyancy frequency. This is the frequency at 
which a fluid particle will oscillate if it is vertically 
disturbed from its equilibrium location. 

The task now is to use the new equations of motion 
to derive an equation that describes the vertical motions of 
the water particles in time and space. This 1s done by 


combining the time derivative of the continuity equation 


(12) and the directional derivatives o£ the horizontal 


momentum equations (9,10). 
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where V7, =0°/dx?+07/dy? is known as the horizontal 
Laplacian operator. 
Equation (14) is still an equation of multiple 
variables. The elimination of pressure from the equation 


can be accomplished by first taking the time derivative of 


equation (11) and using equation (13). 
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Now the horizontal Laplacian of equation (15) is 
taken and substituted into the time and depth derivative of 


equation (14). 
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derivatives of equations (9,10) and the continuity equation 


(12) yields a differential equation for the vertical motion. 
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where V’ is the three dimensional Laplacian 


operator. 


Equation (17) is a second order, homogeneous, and 
separable differential equation. An assumed solution to 


this equation is a periodic wave. 
w= W, meres ( ies: ) 


If the assumed solution is substituted into the 
differential equation and the terms arranged, we get the 
governing equation for internal waves from which the 
dispersion relation is easily shown. 
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Assuming the vertical structure is of the form 
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The dispersion relation gives equation (20) two 
types of solutions. The first type is an exponential 
solution and occurs when (@W>N) or (@<f). In order for 
the equation (20) to have a wave solution, the frequency 
must be bounded by f<@<N. This is the necessary 
condition for freely propagating internal gravity waves to 


exist. 


b) Internal Wave Reflection 

Internal waves are difficult to visualize for most 
people. To help in  theim@gevisualazatieom “Morr (19c1) 
suggests looking at a wave packet with crests and troughs 
along planes normal to the page and inclined with respect to 


the (x,z) axis as shown in Figure 1. The phase velocity is 


— 


in the direction of k (wave number vector) normal to the 
Crests, and the greupevel 6am Caicmeara le togehe sence sece 
The net result is that the wave packet appears to slide 
Sideways. The angle of the wave number vector (k) relative 
to the horizontal can be found from: 
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This results ina steep angle of phase propagation 
when the internal wave frequency is near the inertial 
frequency (w= f); and in a flat angle of phase propagation 
when the internal wave frequency is near the buoyancy 


frequency (w= N). The group velocity is directed 90° from 


the phase velocity with vertical component of opposite sign. 
This means that for near-inertial frequencies, (w= f), the 
packet of energy will propagate at near horizontal; and at 
near-buoyancy frequencies ( =N), the energy packet will 


propagate at near vertical. In both cases, linear theory 


holds that the group velocity will approach zero as the 


an@ 


internal wave frequency approaches either the w=f or the 
@=N limit. 

When internal waves’ reflect, the angle of 
incidence relative to the horizontal equals the angle of 
reflection relative to the horizontal. This is an important 
distinction from other types of waves that preserve their 
angle relative to the reflecting surface. Figure 2 
illustrates the geometry of internal wave reflection. Lae) 


the figure, 2z’ is the vector normal to the reflecting 
surface, Q@ is the slope of the reflecting surface, x' is 
the vector parallel to the reflecting surface, g is the 
gravity vector, and @ is the angle with respect to the 


horizontal at which the internal wave group propagates. 
Because of the horizontal angle reference, 
internal waves have three types of reflection: subcritical, 
Supercritical, and critical. Subcriticalereflection OGemnas 
when the slope of the reflecting surface is less than the 
angle of energy propagation for the wave packet (referenced 
to the horizontal). These waves will be forward reflected 
towards shallower waters. Supercritical reflection occurs 
when the reflecting surface has a slope steeper than the 
angle of wave energy propagation. These waves will be back 
reflected towards deeper water, where they will likely 
dissipate in the open ocean (Munk 1981). Cita] Gat 


reflection occurs when the slope of the reflecting surface 
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1s equal or nearly equal to the angle of wave energy 
propagation. At critical reflection angles, linear wave 
theory develops a singularity (Huthnance 1989). Energy 
density becomes infinite and group velocity falls to zero, 
violating the underlying linear assumptions. The result of 
this may be the formation of non-linear effects such as a 
bore or other hydraulic jump traveling up the slope of the 
reflecting surface (Slinn and Riley 1996). 

2 Non-Linear Internal Wave Theory 

Internal waves are not fully described by linear wave 
theory. Numerous observations have shown that many of the 
features of internal waves can only be explained by non- 
linear theories (Garrett and Munk 1979). Non-linear wave 
theory iS a complicated subject, so only a very. small 
portion of it will be discussed here. For a more complete 
summary of non-linear internal wave theories, see Ostrovsky 
and Stepanyants (1989). 

A soliton is a solitary non-linear wave. A soliton can 
exist either as a surface gravity wave or as an internal 
wave. They have been observed in both forms in nearly every 
part of the world’s oceans (Ostrovsky and Stepanyants 1989). 
A balance between the non-linear steepening of the wave and 
dispersive effects results in the form of a soliton being a 
Single hump propagating at a phase speed that increases with 
the amplitude of the soliton. Various non-linear theories 


can adequately describe the soliton. These include, but are 
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not limited to, the KdV, modified KdV, Joseph, and Smooth 
Stratification Model. Perhaps the most often used theory to 
describe the non-linear internal waves is the Korteweg-de 
Vrres (Ke mecso butaome A soliton is one of two possible 
solutions to the Korteweg-de Vries equation, the other being 
a train of cnoidal waves. However, caution must be used 
when applying these theories to solitons or bores as most 
non-linear theories assume that the non-linearity will be 


weak (Ostrovsky and Stepanyants 1989). The rataemmor 
nonlinear to dispersive effects (a? /H? | is used to determine 


how the form of a soliton will change over time where a is 


the amplitude of the solitary wave, H is the depth of water 
the wave is in, and A is the wavelength of the solitary 


wave. When the ratio has a value of less than 16, the non- 
linear effects can balance the dispersive effects and the 
soliton maintains a constant form, however when the ratio is 
greater than 16 the non-linear effects overpower the 
dispersive effects and lead to a steeping wave front 
eventually causing a hydraulic jump (Kundu 1990). Perhaps 
the simplest model of a hydraulic jump is the 2-layer 


hydraulic jump model used by Holloway (1987). Figure 3 


shows a diagram of this model. In the figure, h, is the 


initial bottom layer thickness with velocity uj, and h, is 


ime 


the bottom layer height after the hydraulic jump with 
ve keecmty um 
Holloway (1987) defines the conditions required for a 


hydraulic jump in terms of the Froude number (F) 


2 
2 u, 


F? =— 
gh, 





(23) 


where g’=(p, —p,)g/p, 
He defines a critical F where it becomes possible for a 
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with h,H,u,,u, from Figure 3, which shows a downward jump in 


the upper layer if the front is considered stationary. For 
a moving jump, flows relative to the jump speed should be 
used in equations (23) and (24). He found that, hyememiic 
jumps occurred approximately three hours after F became 
emreneal auvay value of 2 much) larger thane) Generaesmg a 
larger hydraulic jump than would have been possible if F had 
just barely been super-critical (Holloway 1987). The 
largest downward jumps were found to exist when there was a 
thin upper layer present over a weak lower layer flow. 
Surface jumps are common throughout the world, from 


Sspillways on dams to jumps near large rocks in quickly 
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flowing rivers. A bore 1s a propagating hydraulic jump. 
These are less common than stationary jumps. Surface bores 
have been observed in a few rivers and inlets while internal 
bores have only rarely been observed (Halpern 1971, Holloway 
moe7) Tidal bores in rivers are formed by interactions 
between the tidal wave and the narrowing river channel. 
Basically the tidal wave will grow as the tidal energy is 
confined to a smaller and smaller area until it becomes 
unstable and develops into a hydraulic jump (tidal bore). 
These tidal bores tend to be largest during the spring tides 
(Pugh 1987). The tidal bore on the Severn in England is an 
excellent example of a surface bore. 
ae OBSERVATIONS OF NON-LINEAR INTERNAL WAVES ON 
CONTINENTAL SHELVES AND SEAMOUNTS 

Lee and Beardsley (1974) used data collected from 
Massachusetts Bay by Halpern (1971) to compare with a non- 
linear KdV-type wave in an effort to determine the relative 
importance of non-linearity and dispersion on the generation 
of large amplitude internal waves in “the breakdown of 
internal fronts. They found that their model adequately 
described the physical features found by Halpern in his 
observations of tidally generated internal waves. 

Howell and Brown (1985) found soliton packets on the 
Northern California continental shelf during the Coastal 
Ocean Dynamics Experiment in 1981. They concluded that the 


internal solitary waves are generated on the continental 


i 


shelf and evolve into soliton packets as they propagate 
shoreward. 

Holloway et al. (1997) use a numerical solution to the 
KdV equations to model the evolution of an initially 
Sinusoidal long internal wave to represent an internal tide. 
Their model shows a transformation of the initial waveform 
into a series of shocks and solitons as it propagates over 
the continental shelf. Their model was run uSing 
hydrographic data from observations made over the Australian 
Northwest Shelf. 

Stanton and Ostrovsky (1998) reported strong solitons 
off the Oregon coast. Extremely strong solitary internal 
wave packets were observed over a three-week period on a 
very shallow and strongly stratified pycnocline off northern 
Oregon. During periods of the strongest tidal forcing, 
solitons were consistently observed on the leading edge of a 
semidiurnal tidal bore, with pycnocline displacements of up 
to 25 m downward from a 7-m initial depth. The extreme non- 
linearity of these observations are beyond the limits of the 
weak non-linearity assumed by the KdV model. A higher order 
model, the "CombKdV model", is proposed to describe the 
observations. The higher order model, despite its 
assumption of still weak non-linearity, closely matches the 
observed data. 

Eriksen (1998) observed internal wave S[HoRe EMEEOS as 


they reflected off the Fieberling Guyot seamount in the 
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Néaeheeeacl £1 C..0cean: His observations were quite long, 
lasting approximately one year. He observed that the 
records substantially differed from the predictions of 
linear wave theory. Some interesting observations were the 
size of the reflections and the presence of a critical slope 
reflection. The effects of internal wave reflections in the 
primary mode of the EOF's were detectable at 750 m above the 
seafloor. These reflections had a phase structure that 
conformed to the linear theory predictions, but the 
amplitude decayed in proportion with the EOF wave number. 
Intensifications of the wave spectrum with respect to the 
GM79 spectrum were seen at the local critical frequency. 
The near bottom water was well mixed due to large and 
frequent density overturnings caused by the wave 


reflections. 


D. INTERNAL WAVE WORK IN MONTEREY AND OTHER CANYONS 

The Monterey Canyon is located nearly in the middle of 
Monterey Bay, California. Its upper reaches approach the 
coast to within a few hundred meters at Moss' Landing, 
SarerrOriita, making 1t very attractive Beng deep-sea 
experiments and observations. Both Moss Landing Marine 
Laboratories and the Monterey Bay Aquarium Research 
Institute (MBARI) are located near the canyon head and make 
frequent, mainly biological and geological, measurements in 


the canyon. The surface tide in Monterey Bay is of the 


ly, 


mixed semi-diurnal type and essentially co-oscillates 
(Petruncio et al. 1998), with a strong internal tide in the 
Monterey Canyon (Petruncio et al. 1998). 

Broenkow and McKain (1972) observed tidal oscillations 
over a 25-h period from two locations 3 km apart. Their 
measurements revealed semidiurnal isopycnal oscillations 
with maximum ranges of 80 m propagating shoreward at 20 to 
30 "eni/ Ss" 

Shepard et al. (1979) found current oscillations with 
periods of eight to nine hours that were dominated by the 
overtides rather than the semi-diurnal internal tide from 
three 4-day records of current measurements recorded at 3 m 
and, 30 misover the camyon floer. They demonstrate some of 
the difficulties with observations in Monterey Canyon. 
Unlike many of the other canyons discussed by them, the 
Monterey Canyon does not generally have water flow that is 
aligned with the canyon axis. They speculated that there 
may be many more tributary canyons or that the canyon 
topography is just much more complicated than believed at 
the time. As the hydrographic survey conducted by MBARI in 
1998 showed, their thoughts on the canyon topography were 
correct; Monterey Canyon has a very complicated structure. 

Shea and Broenkow (1982) used data collected by 
Broenkow and McKain (1972) and data collected in 1978 and 
1979 to further analyze the internal tidal signal looking 


for evidence of nutrient enrichment. They found an internal 


ie 


fide with herehts ranging “rom 50 to 120e@m and the presence 
of a breaking internal bore 25 m above the bottom near the 
canyon head. The bore was identified using a thermistor 
just offshore of the Moss Landing harbor entrance. 

The Monterey Bay National Marine Sanctuary site 
characterization, physical oceanography section (Broenkow) 
reports internal tidal displacements with magnitudes of up 
to 120 m just 5 km from the canyon head. Other measurements 
are cited from Heard (1992) with 60 m internal tidal 
amplitudes in 185 m of water near the canyon head. 

Petruncio et al. (1998) analyzed data from ITEX1 and 
ITEX2 experiments conducted in Monterey Bay during April and 
October, 1994, respectively. Vertical profiles ohn 
temperature, salinity, and velocity (excluding the bottom 20 
to 30% of the water column) were made at one to two hour 
intervals over 25-h periods during each experiment. They 
found semidiurnal currents in the 20 cm/s range, an order of 
magnitude larger than the estimated barOeropl Gat aeal 
currents. They proposed a link between surface and internal 
tides as they relate to the stratification of the water. 
Shoreward energy propagation was observed during ITEX1 anda 
Seanding wwave (interne lee seucne ) securing BelrE ws. They 
speculate that this was due to differences in the 
stratification present in Monterey Bay during April and 


October when the two experiments were carried out. 


des, 





Their analysis of the internal tide was limited by 
design of the experiments to describing the linear internal 
tide. Based on the assumption that the linear internal tide 
can not negotiate the sharp bends in the canyon, two 
possible generation sites, Smooth Ridge and Steep Ridge, for 
the locally observed internal tidal signal were identified. 

Lafuente (1999) examined the La Linea Canyon in the 
western Alboran Sea. He focused his measurements on the 
shallower areas of the canyon (< 100 m), but did make some 
measurements in the deeper areas (~300 m). He found that 
the deeper tidal currents exhibited large amplitudes (~50 
cm/s) and were phase locked with the more barotropic tidal 
currents near the surface. 

Matsuyama et al. (1993) measured the near-bottom 
currents in the Suruga Trough in Japan. They found the 
currents were dominated by the M2  semi-diurnal tidal 
constituent, which was five times the strength of the next 


strongest constituent. The near-bottom currents seemed to 


be phase locked near 0° with the sea level while the 


barotropic tidal currents were 90° out of phase with the sea 


level. 

Takeuchi and Hibiya (1997) made observations near the 
Suruga Trough in Uchiura Bay. They identified a generation 
Site for the internal tide in Uchiura Bay and conducted 


measurements in the area. They discovered that the 


ay 


Character of the internal tide in the bay was dependent on 


the non-linear dynamics at the generation site. 


ou 
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II. DATA COLLECTION 


A. MOORED 

Data were collected in the Monterey Canyon over a 34- 
day period from 05 August 1997 until 10 September 1997. ‘Two 
moorings, BB and NB, were deployed approximately 2.8 km 
apart from each other after a brief bathymetric survey was 
conducted (Figure 4). The moorings were lowered to within 
10 m of the bottom before being released. Both were 
intended to be close to the middle of the canyon axis, but 
high resolution bathymetry obtained the following summer by 
MBARI revealed that BB was on the axis at the apex of a 
curve in the canyon, while NB was off the principal axis on 
the southern flank of the canyon (Figure 5). MBARI obtained 
their high-resolution bottom topography using a Simrad 300 
bottom mapping system. It has a vertical resolution of 0.1% 
of the water depth and a horizontal resolution of 2% of the 
water depth. 

Both moorings were short, taut moorings consisting of a 
railroad wheel anchor, two acoustic releases and an upward 
looking, self-contained ADCP deployed in a foam buoy. An 
Alpha Omega thermistor was attached to the metal frame 
protruding from below the buoy. 

One mooring, NB, supported a RD Instruments 300 kHz 


narrowband ADCP ina 48" diameter buoy. The ADCP transducer 


Ze 


was 6 m above the seafloor in ~390 m of water and resolved 
currents in 4-m bins with the deepest bin centered 12 m 
above the bottom. Velocity data were calculated by vector 
averaging the results from 270 pings spaced evenly over 
3-min intervals. RDI reports a nominal standard deviation 
Gue to instrument error of 0.8 cm/s for a 270-ping average 


for the narrowband ADCP. The Alpha Omega temperature 
recorders employ a YSI 46006 thermistor with +0.03 °C 


stability. The instruments were calibrated by the Oregon 
State University mooring group immediately prior to the 
experiment. The clocks maintained accuracy to within 40 s 
over the course of the experiment. The thermistor was 
located 6 m above the bottom and recorded instantaneous 
temperatures at 1-min intervals. 

The second mooring, BB, supported a RD Instruments 300 
kHz broadband Workhorse ADCP ina 37" foam sphere. The ADCP 
was 6 m above the bottom in ~333 to 337 m of water. The 
depth of the BB mooring was determined by comparing the 
pressure sensor in the ADCP, the fathometer readings taken 
from the R/V Pt. Sur during deployment, and DGPS- determined 
location relative to the high-resolution bathymetry data 
collected by MBARI in the spring of 1998. For the purposes 
of this analysis, the depth of the BB mooring anchor will be 
taken as 333° mM: It resolved currents in 4-m bins with the 


deepest bin centered 12 m above the bottom. Velocity data 
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were determined by vector averages of 12 pings spaced evenly 
Gver ies Minwantervals: RDI reports a nominal standard 
deviation due to instrument error of 0.8 cm/s for a 12-ping 
average for the Workhorse ADCP. The BB thermistor recorded 
an instantaneous temperature at 6 m above the bottom at 1- 
min intervals. Both ADCP’s recorded velocity data in earth 
coordinates and an instantaneous value of pitch and roll in 
each ensemble header at 1.5-min(BB) or 3 min(NB) intervals. 
The BB and NB moorings were placed in such close 
proximity, and used similar sampling schemes, to provide 
redundancy in case of failure of one of the moorings. As 
will be seen, both instruments worked well and the records 
show a high degree of spatial variability in the canyon. 
ADCP’s use backscattered sound to measure velocity. 
The Doppler effect is used by transmitting sound at fixed 
frequencies and listening to the returning echoes from sound 
scatterers in the water. The sound is scattered from the 
small plankton and particles in the water column, which are 
assumed to have, on average, the same horizontal velocity as 
the water that they are in. The frequency reflected by the 
sound scatterer is shifted lower when the scatterer is 
moving away from the ADCP and shifted higher when the 
scatterer iS moving towards the ADCP. As the Doppler shift 
is only sensitive to along-beam motion, the components of 


flow along several radial directions, and the orientation of 


Zo 


the transducers, must be measured to estimate the full 
velocity vector. 

The ADCP’s used in this experiment have four acoustic 
beams which can determine three orthogonal components of 
velocity, with redundancy. The four transducers are 
arranged in a Janus configuration that allows beams on 
opposite sides of the square to act aS a pair. The method 
assumes the velocity field is horizontally homogeneous over 
the area ensonified by the four transducers. Due to the 
large influence of mooring motion on vertical velocity 
estimates, only horizontal components of velocity were 
resolved in these data. 

We considered that the mooring motion may have 
introduced possible errors in the horizontal velocity data 


as well. While the maximum pitch and roll recorded by the 


tilt sensors "reeords onwene two ADGE "Ss was 325 > welle@wiehin 
the + 15° operating range specified by RDI, the maximum 


angular acceleration was 2.3° per min for the BB ADCP. 


Knowing that the large angular accelerations experienced by 
surface moorings can contaminate ADCP velocity measurements, 
several messages were exchanged with RDI on the response of 
the tilt correction algorithms to rapid (on the order of 
minutes) tilt angle changes. At this time we do not have 
full characterization of the reaction time of the tilt 


correction. 
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Velocity data degraded by erroneous tilt values are 
typically seen as rapid accelerations/decelerations that are 
coherent throughout the measured water column. Even though 
no coherent changes were present in the records during the 
strong tilting events, two further steps were taken to 
convince us of the validity of the measured velocities. To 
get an order of magnitude estimate for the possible errors 
in the data due to any uncorrected tilt angles, calculations 


were performed on a typical tilt excursion. Using a maximum 
tilt angle of 3° in a 1.5-min interval, the maximum speed 


deviation found to exist is less than 7 cm/s. An analysis 
of the mooring motion confirmed that the measured tilts were 
consistent with the predicted instrument deflections based 


on the recorded currents. 


B. SHIPBOARD 

Four 12-h CTD time series were taken from the R/V Pt 
Sur. The length of each series was limited by the need to 
turn the ship over to the operators of: the opposing 12-h 
watch. CTD stations are marked on Figure 4 as square boxes, 
the station between the two moorings is known as "Shallow" 
and the other station is known as "Deep". Measurements were 
taken at Shallow on 09 and 21 August 1997 and at Deep on 10 
and 20 August 1997 using a 9/11 SeaBird CTD along with a 
SeaTech light transmissometer with a 25 cm path length. CTD 


casts were conducted at Shallow approximately every 30-min 


au 


from 0708Z - 1848Z 09 August 1997 and again from 1029Z - 
2140Z. 21 August 1997. 

Current data were recorded by a 150-kHz ADCP (VM-ADCP) 
mounted in the hull of the R/V Pt Sur while CTD casts were 
being conducted. The VM-ADCP calculated average currents 
over either a 3-min or a 1l-min interval. The ship's heading 
and velocity, obtained from gyrocompass and differential GPS 
inputs respectively, are needed to obtain earth-referenced 
water velocities. Some of the VM-ADCP data was found to be 
unusable due to problems with either the GPS and/or 


gyrocompass interface and were not used in this analysis. 


Cs TIDE GAUGE 

Sea level data were obtained from tide gauge No. 941- 
3450, located in Monterey, California and maintained by 
NOAA . The gauge can be found in the harbor maintenance 
Supervisor's office at the seaward end of Municipal Wharf 
No. 2 in Monterey. 

The sea level data are collected using the Next 
Generation Water Level Measurement System which uses an 
acoustic measurement device rather than the older float type 
system. The water level is measured by sending a shock wave 
of acoustic energy down a half-inch diameter pipe and 
measuring the travel time using a calibrated reference. 
Another change from older stilling well float measurements 


1s in the well itself. The new protective well is more 
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exposed than the stilling wells were, allowing more of the 
wind waves and chop to affect the sea level measurements. 
Acoustic measurements are recorded every 6 min from the 
average of a set of 181 1-s interval measurements. Data 
recorded are the mean sea level, standard deviation, and the 
number of outliers with a resolution of 0.01 feet. Hourly 
data, obtained from NOAA'S web site, were used in all the 
figures presented here, with the exception of Figure 17 and 
Figure 18 which used the times of high and low sea level to 


the nearest 6 min. 


D. ANCILLARY 

Various other data have been used to corroborate 
assorted findings. These include: AVHRR data from NOAA 
satellites, data from various buoys in Monterey Bay placed 
by MBARI, and surface currents derived from a CODAR type HF 


radar system operational in Monterey Bay. 
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IIIT. ANALYSIS AND RESULTS FOR FULL TIME SERIES 


A. PRINCIPAL AXES 

Data from the two ADCP's was recorded in earth 
referenced coordinates, but the currents from the two 
moorings do not lie along either the North/South or the 


East/West axis. To further complicate data analysis, the 
up- and down-canyon currents are not aligned 180° apart from 


each other. 

Two methods were used to find the "best" way to 
calculate the principal axes for the current data. The 
first method calculates the principal axes by assuming the 
up-canyon direction is in the NE quadrant and the down- 
canyon direction iS somewhere in the other three quadrants. 
Once the principal axes are calculated the current data are 
projected onto the derived axes (Figure 6). Data in the NE 
and SW quadrants are projected onto the positive and 
negative A' axis respectively, while data in the NW and SE 
quadrants are projected onto the negative and positive A" 
axis respectively. The second method created two least- 
squares-fit lines from current data in the u>0 and u<0 
regions separately. The results of the two principal axes 
calculations are show overlain on the current data from BB 
bin #1 (Figure 7) and a comparison between the original 


speeds and the speeds along the new along-canyon axes was 
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Completed (Figtuwems | Negative values indicate westward 
flow in the original data and down-canyon flow in the 
rotated data. It is clear that the first method is superior 
to the second, replicating virtually all of the variance in 
speed seen in the original data. 

Figure 9 and Figure 10 show velocity scatter plots for 
four depth bins at NB and BB respectively with the derived 
principal axes from method one shown. The calculated near- 
bottom along-canyon axes at both mooring locations are in 
general agreement with the expected directions given the 
topography at the mooring sites (Figure 5). The principal 
axes directions at the BB location show more variability 
with depth than at the NB location; this is attributed to 
the location of the BB mooring near a sharp bend in the 
canyon topography. The component of flow along the 
principal axis in the NE quadrant is subsequently referred 
to as "“up-canyon", while the component of flow along the 
principal axis in the NW or SW quadrant, depending on bin 


number, is referred to as “down-canyon". 


B. POWER SPECTRA 

Spectral analysis was performed on the current records 
obtained from the BB ADCP. The spectra were calculated by 
averaging the fast Fourier transform results from 9600-point 
Hamming windowed pieces (corresponding to a 240 h time 


period) overlapped by 50%. In addition to the piece 
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averaging, each bin's spectrum has been frequency averaged 
from Q.3-cph onward with a linearly increasing average 
length. As expected, there was a strong spectral peak at 
the M2 frequency (12.42 h) (Petruncio et al. 1998). Figure 
11 shows the kinetic energy spectra for all 29 depth bins 
from the mooring, with each bin separated on the y-axis by 
one decade. In addition to the prominent semidiurnal peaks, 
numerous overtides are clearly seen, including peaks at the 
M3, M4, 2MK5, M6, and 3MK7 frequencies. We have examined 
the spectra all the way out to the Nyquist frequency and 
there are no discernable humps or change in slope near the 
Brunt-Vaisaéla frequency. 

The kinetic energy spectra for the top (bin #29) and 
bottom (bin #1) depth bins, 124 and 12 m above the bottom, 
respectively, were compared to the kinetic energy spectra 
for open ocean internal waves developed by Garrett and Munk 
(GM79) (Figure 12). The kinetic energy for both depths is 
higher than predicted by GM79 and the near-bottom kinetic 
energy is higher than the shallower kinetic energy. The 
near-bottom spectrum's peaks at the M3, 2MK5, and M6 
frequencies are nearly the same level as the peak at the M4 


frequency for the shallower spectra. 


Oe TIDAL ANALYSIS 
Tidal analysis, resolving 35 different tidal 


constituents, was performed on the current records. It was 
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found that ae kinetic energy of the deeper currents was 
spread out more evenly over the various constituents than 
the shallower currents. 

Table 1 summarizes the differences between the depth 


bins for selected tidal constituents. The strength of the 


odd numbered shallow-water constituents in many bays and 


narrow estuaries are typically not very strong in semi- 


diurnal tidal regimes, however mixed tidal regimes may have 


the fifth-diurnal tides enhanced (Pugh 1987). In the 


canyon, which may be thought of as like a river, the £1ifth- 


diurnal tide is enhanced, stronger than the M6 at some 


depths. The third and seventh-diurnal tides are also 


enhanced in our data. Curiously the third-diurnal tide has 
a minimum in the middle depth bins while the other odd 


numbered constituents have a maximum in the middle depth 


bins. 
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D. TIME SERIES 


Figure 13 shows the thermistor records from the Alpha 


Omega thermistors on the NB 


with the sea level from the Monterey Bay tidal 


(green) 
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and BB 


(blue) 


moorings 


station 


overlain. Initial inspection of Figure 13 immediately 
reveals the presence of periodic rapid temperature drops at 
both the NB and BB locations, separated by more gradual 
periods of warming. It appears that the temperature drops 
are due to the passing of an internal tidal bore from 
farther down-canyon, and will be discussed further. 

The temperature drops at the two mooring locations have 
different characters (Figure 14). As seen on O09 August 
1997, the BB records clearly shows two large temperature 
drops, whereas the NB record shows four smaller temperature 
drops. The difference in the two records may be due to the 
locations: of each mooring relative to the canyon topography. 
The MBARI bathymetry data indicates that the NB mooring was 
south of the canyon axis on a bank ~40 to 50 m above the 
floor of the canyon axis. The BB mooring appears to have 
been just off the center of the canyon axis. The actual 
locations of the two moorings made the BB mooring a natural 


selection for a more in-depth analysis. 
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IV. DESCRIPTION OF THE TIDAL BORE 

The most prominent signal in the BB mooring temperature 
ane SLOci ty records swadesbe wre Lemme witOna. Dr10Orl aS thoteot 
a tidal bore. An accurate and reliable method to identify 
the passing of a bore was sought. It was found that the 
tilt values recorded by the ADCP pitch and roll sensors 
showed sudden changes that were highly correlated with the 
rapid temperature drops (Figure 15). It was also found that 
an increase in the up-canyon current accompanied the large 
tilt excursions and temperature drops, as shown for the 
pitch and near-bottom currents on 09 August 1997 (Figure 
16). This figure shows the typical relationship among these 
three variables. 

It seems unlikely that the mooring motion caused by the 
bore passage seriously contaminated the measured velocities. 
Any velocity oscillations at the period of the mooring 
motion are small (Figure 16) and there is no coherent high 
frequency signal throughout the water column that matches 


the period of the experienced tilts (Figure 21). 


A. TIMING RELATIVE TO SEA LEVEL 

The bore signal changes throughout the 34-day record 
(Figure 13). To investigate a possible relationship with 
the barotropic tidal signal, the timing of the bores 


relative to sea level as measured by the Monterey tide gauge 


on 


was examined. For the entire 34-d record, the onset of each 
bore was visually identified in the temperature record and 
the time difference between events was calculated. The time 
difference varied between 1 and 24 hours but in general fell 
between 10 and 14 hours with a mean of 11.9 h, slightly less 
than Sthe “pericdmwort the M2 “constrluent 112242. Hh): The few 
time intervals exceeding the K1 period and the few extremely 
short intervals may result from errors in the visual reading 
of the temperature records. 

During 01 to 09 September 1997, large signals in 
instrument tilt, near-bottom temperature, and near-bottom 
velocity occurred with great regularity. This time period 
was examined in more detail. The times of bore passages 


were estimated from the tilt records as the time when the 
angle of pitch exceeded + 1°. Over these nine days, the 


bores had an average separation interval of 12.3 h, with a 
maximum of 13.4 h, a minimum of 11.3 h, and a standard 
Geviatten Of 0.6 n= Over the same time period, the high 
tide had a mean interval of 12.3 h, a maximum of 13.4 h, a 
minimum of 11.4 h, and a standard deviation of 0.6 h (Figure 
17). Over these nine days, the time interval between a bore 
passage at mooring BB and the previous high tide had a mean 
of 8.6 by with a maximumnm@em 9.5 h, a mmnimum of 7.8 Hh, anda 


Standard deviation of 0.6 h (Figure 18). 


SS 


To examine the hypothesis that the temporal variation 
in the strength of the internal bores may be related to the 
spring/neap cycle, the ADCP pitch was used as a proxy for 
identifying the passage of the bores. The pitch record is 
compared to the tidal range, calculated as the range from 
either high tide to low tide or low tide to high tide, over 
the entire 34-day record (Figure 19). It is immediately 
obvious that two of the three sets of strong bores (~8/6- 
8/11 and ~9/1-9/9) occur when the tidal ranges are close to 
their minimum while only four strong bore signals occur 
during the three-day period when the tidal ranges are 
greatest (~8/16-8/19). From 21 to 24 August, when the tidal 
ranges are again near aminimum, only two strong bores are 


noted. 


B. VERTICAL STRUCTURE 

Since data from ~12-h time-series of CTD casts were 
available at the Shallow CTD station on 09 and 21 August 
1997, these two days will be examined in more detail. These 
two observation periods seem to fall in different regimes 
for the bore. While a significant bore passed through this 
region of the canyon on 09 August, only weak, or no, bores 
were seen on 21 August. 

ane 09 August 1997 

Figure 20 shows all of the density profiles derived 


from the CTD casts conducted at ~30-min intervals at the 
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ShalloweGEiDn stabmen from 07 0sZeunti wie 4 ze in Saddatsien , 


the average density profile and the N profile calculated 
from it (using Az of 1m) are also shown. Perhaps the most 


notable feature in this figure is the sharp pycnocline. At 
~200 to 250 m depth, near the depth of the upper bins from 
the moored ADCP, the temporal variance is much larger than 
near the bottom. 

Figure 21 shows the BB and vessel-mounted ADCP total 
speed data (both 1 min and 3 min), isopycnal depths derived 
from the CTD data, and the temperature data from the moored 
Alpha Omega thermistor during the time period of the CTD 
station. This figure shows that CTD casts began shortly 
abeen ele Pao slilgeol va WOre ane ene SMOmely erOme tellers tema 
of the next bore. The internal tide signal is clearly seen 
in both the moored ADCP records and the vessel mounted ADCP 
records, showing a reasonable match at the interface between 
the two measurements. Figure 22 is a plot of the east and 
north currents from the BB ADCP at a depth of 217 m and from 
the vessel mounted ADCP at a depth of 215 m. The two ADCP's 
show better correlation in the phase than in the magnitude 
of the currents. This difference is attributed to the 
spatial offset of over 1-km between the BB ADCP and the 
location of the CTD casts. Even though the vessel mounted 
ADCP was not directly over the BB mooring, it was mainly 


over the canyon axis farther down canyon. Given the 
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vessel's location, if the bore affected the upper water 
column, some evidence of an offset bore passage should be 
seen in the vessel mounted ADCP record. As seen in Figure 
21, there 1s no evidence for bore effects above 209 m. 

The first bore on this day, seen at ~0700, is shown in 
detail in Figure 24. It is easier to see in this figure how 
the CTD casts begin shortly after the commencement of the 
07-00 = boner This particular bore is associated with one of 
the largest temperature drops recorded during the entire 34- 


day record. The initial temperature drop for this bore is 


“ies )C am 1 min: The maximum range for this bore is ~- 


2 eT) ee ae ihe indtlial Grop 15.0402) pone ener teaman 


temperature decrease and occurS over just one sampling 
interval of the thermistor. Maximum current speeds close to 
the bottom are up-canyon in the 45 to 55 cm/s range and 
occur within 30 min of initial temperature drop. Currents 
accelerated in the up-canyon direction with a maximum value, 
~0.24 cm/s’, at 12 m above the bottom coincident with the 
time of the initial temperature drops. Starting at ~0730, 
thin "finger-like" regions of increased current speed 
extending upward 80 to 90 m from the bottom are present for 
~20 min. A maximum recorded current of 69 cm/s occurred ~40 
m above the bottom (bin #8, 32 m above the transducer) 
approximately 35 min after the initial passing of the bore. 


These high velocity events in time and depth are associated 
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with current relaxations between the short current pulses 
higher in the water colum. Figure 23 shows that the 
current actually relaxes rather than totally reversing 
between each of the high velocity events. 

The second bore on O09 August occurred at ~1845 and has 
Similar velocity structure as the first (Figure 25). The 


temperature signal for this event is not as clear. The 


maximum temperature drop associated with this bore is 1° C 


in 50 min with a maximum single step of just 0.25° C 


occurring near the end of the temperature drop. This single 
step represents 25% of the total temperature range. The 
maximum current speeds near the bottom are in the same range 
as the earlier bore (45 to 55 cm/s) and occur with the same 
EtmMing (within 30 min of initial temperature drop). 
However, this bore lacks the high-speed current pulse seen 
in the earlier bore at 40 m above the bottom. The maximum 
acceleration seen at 12™m above the bottom 1s  —U]l2 cm/s 
approximately 60 min after initial bore passage, and does 
NOeProcecumrae = time of initial bore passage which has an 
acceleration of ~0.1 cm/s’. Tie finger-like” =struetures 
above this bore are more uniform in time and space and they 
are easier to identify than in the earlier bore. 

In addition to the Signal associated with the two 
temperature drops, a series of pulses of high speed flow 


(~50 cm/s) occurs ~60 m above the bottom between 250 m and 
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270 m depth from ~1400 until ~1800° (Figure 6)— These 
strong pulses enhance the generally down-canyon flow 
G¢eurring during this time periods The pulses oeceurwwith an 


average period of 2.5 min and gradually increase in strength 


over 60 min before dying away. The period of these pulses 
1s longer than the period of the buoyancy oscillations (~1 
min) at these depths, making it possible for these 


oscillations to be freely propagating internal waves, 
perhaps generated from a small ridge or other topographic 
feature up-canyon of the BB mooring site and advected down- 
canyon by the lower frequency background flow. 

Bae 21 August 1997 

Figure 27 shows the density profiles calculated from 
the CTD casts conducted at the Shallow CTD station from 1029 
Wimied ee) tem 62 l  AUGUSEN lo 7. The difference in density 
structure between 21 August (Figure 27) and 09 August 
(Figure 20) is immediately apparent. The water column is 
more continuously stratified on 21 August and the pycnocline 
much less distinguishable. In fact, comparison with the 
data collected at the deep CTD station indicates that much 
of the change took place between 09 and 10 August. Buoy 
measurements indicated that winds reversed from equatorward 
to poleward on O07 or O08 August (depending on location), and 
HF radar-derived surface currents across the mouth of 
Monterey Bay show a Similar reversal on 08 August. While 


cloud cover limited availability of AVHRR data between 08 
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and 11 August, the 07 August image shows the characteristic 
upwelling pattern, whereas 12 August 1s uniformly warm. 
Winds are again equatorward, although not particularly 
SEONG) during 2 no” SaarAuGguseE, followed by a weak 
relaxation during 16 to 18 August. Surface temperatures are 
high throughout this period. Previous work (Rosenfeld et 
al. 1994, Ramp et al. 1997) has shown that a warm meander, 
of significant vertical extent, moves shoreward into the bay 
during relaxations of the upwelling-favorable winds typical 
FOr me this came Geof — year. The observed changes in 
stratification could well be due to such a mesoscale event. 
EFarlier observations of this phenomena would indicate that 
the density adjustments take place within one day of the 
changes in the winds, whereas an interval of at least two 
days must have elapsed if this mechanism is responsible for 
the changes observed in these data. 

Figure 28 is a contour plot of the moored ADCP data 
from 21 August. Unfortunately, the shipboard ADCP was 
experiencing difficulties on that day making full water 
column measurements impossible. However, since on O09 August 
when a particularly strong bore was noted, the effects did 
not extend above 209 m, it is reasonable to assume that not 
much of the signal of interest is missed by not having VM- 
ADCP data on 21 August when the bore appears much weaker. 
On this day CTD casts began between bores and continued 


through the passing of a bore at ~1345. CTD casts continued 


4 4 


well after the bore's passage and ended between bores with 
the velocity and density structure in much the same state as 
when they began. Figure 30 shows the velocity vectors from 
the BB ADCP. 


The temperature drop associated with this bore is very 
weak, taking just over 2.5 h to fall 1° C with the maximum 
SkuC eosgunime OL: GJtiStemonwO> er im slo =min (Se) Ofmmeocral 
temperature drop). The maximum current speed 12 mab is 44 
cm/s, aS compared to 55 cm/s in the 0700 bore of 09 August 
(Figure 29). It occurs ~45 min after the temperature starts 
to drop. The maximum acceleration at 12 m above the bottom 
is ~0.05 cm/s’ and does not occur at the time of initial bore 
passage which has an acceleration of ~0.03 cm/s’. These 
accelerations are barely distinguishable from the background 
flow, which has accelerations from 0.01 to 0.03 cm/s’. The 
ADCP tilt records show little activity in response to this 
bore (Figure 15). Although the tidal bores on 21 August are 
Significantly weaker than on 09 August, their vertical 


Sige tceers Sam) lars 
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Vv. DISCUSSION 

The kinetic energy in the internal wave band in the 
Monterey Canyon is widely distributed among tidal 
constituents. These include several of the shallow water 
constituents that are a result of the non-linear 
interactions of the internal tide with the ocean bottom 
(Werner 1992). As shown in Figure 11, the energy is more 
concentrated at the M2 frequency higher in the water column 
and spreads out more into the overtides as the bottom is 
approached. Sharp temperature fronts and strongly bottom- 
intensified velocity signals occur approximately at 12-h 
intervals. The time interval between bores, averaging 11.9 
h over the whole 34-d period and 12.3 h during 01 to 09 
September, does not exactly match that of the M2 constituent 
(12.42 h) nor does it maintain a constant relationship to 
the maxima in sea level, indicating that a mix of factors 
comes into play. It 1S not clear where and how these 
features form, why they vary in strength over the 34-d 
record, and what dynamics control them. 

We have described the bottom-intensified signal as an 
up-canyon propagating internal bore rather than as a soliton 
traci: These two seemingly different non-linear phenomena 
can be thought of as a continuum (Munk 1981, Ostrovsky and 
Stepanyants 1989), an initial disturbance that becomes 


bottom-intensified through critical reflection (Thorpe 


4‘7 


1992), steepening of the wave front into a bore, and finally 
evolving into a soliton train (Ostrovsky and Stepanyants 
1989) that breaks in the shallow waters near the canyon 
head. To identify where in the continuum the events 
described in this thesis lay, the characteristics of each 
stage in the evolution were compared to our observations. 
Evidence to support the claim that the events described 
here have fully evolved into a soliton train could not be 
Pee! . First, solitons have usually been observed on a 
density interface near the surface or mid-water column 
(Ostrovsky and Stepanyants 1989, Stanton and Ostrovsky 
1998), and on occasion near the bottom in highly stratified 
water (Bogucki et al. 1997, Holloway 1987). The events 
described here occur near the bottom in weakly stratified 
flow. Second, fully developed soliton trains exhibit a rank 
ordered structure (Ostrovsky and Stepanyants 1989) that is 
not clearly present in our data. Third, the stronger events 
described in this thesis are associated with strong up- 
canyon current accelerations and temperature drops that 
indicate a frontal passage, leaving the fluid in an altered 
State, rather than the temporary changes wrought by a 
soliton. ~ For these reasons and the analogy with tidal bores 
on rivers, these events should be regarded as internal bores 
rather than as a train of solitons. Ultimately, bores and 
solitons are closely linked and can be thought of as steps 


it) the evolution ofr linear intermal waves co turboulence. 
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The finger-like projections of high velocity current 
extending above the bottom layer could be turbulence, high 
frequency internal waves, or even solitons associated with 
the bore front, but without downstream measurements to 
record their evolution and higher temporal resolution 


measurments this remains unresolved. 


A. BORE FORMATION 

Two types of scenarios leading to the generation of 
internal bores have been proposed. One, inspired by 
Halpern's (1971) observations in Massachusetts Bay, consists 
of three phases: formation of an internal front due to 
blocking of the stratified flow by a topographic feature, 
non-linear steepening of the front, and the formation of new 
waves due to the interaction of dispersion and non-linearity 
(Lee and Beardsley 1974). The second involves steepening of 
a large amplitude internal tide as it shoals (Holloway 
ro) Either, or both, of these mechanisms could be at 
work in generating the bores described here. 

Petruncio et al. (1998) have previously shown the 
presence of a large amplitude internal tide in Monterey 
Canyon due to energy concentration by the near-critical 
sloping bot tom and the narrowing canyon walls. 
Furthermore, based on their observations made well above the 
bottom boundary, which compared well with linear theory, 


they suggested that the Monterey Canyon internal tide is 
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generated outside the confines of Monterey Bay. They 
proposed two possible generation sites based on 
characteristic ray tracingmof™the M2 constituent ikal 
identifying these sites, they assumed that while the 
internal tide may be steered by the canyon, the signal would 
not propagate through the numerous sharp bends and curves in 
the canyon bathymetry without significant modification, and 
so the internal wave they observed was likely generated ina 
wedge defined by the north and south wall orientation at the 
canyon head. Data presented in this thesis, from an area 
near a tight bend in the canyon, show that strong non- 
linearities of the internal tide may be associated with 
small-scale topographic features. 

Holloway (1987) examined a series of bores propagating 
onto the Australian North West shelf, and using the 
assumption of two-layer flow, was able to describe the 
conditions required for their formation using Froude 
numbers. The Froude number for a particular layer can be 
calculated using equation (23) and the critical Froude 
number for a bore to form can be calculated using equation 
(24) with appropriate adjustments for the lower versus upper 
layer. The situation observed near the bottom of Monterey 
Canyon does not easily lend itself to analysis in terms of 
layers, nor is it possible with measurements at one fixed 
location to determine the speed of the bore since the 


measured velocities are a combination of the flow relative 
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Gemethe hydraulic jump andthe motion offethe jump ittscit: 
While theoretically over a flat bottom with layered flow, a 
bore could be moving in either direction, while still 
registering up-canyon flow coincident with temperature 
drops, in stratified flow over the sloping bottom an up- 
canyon-propagating bore seems much more likely (Figure 31). 
Conceptually, the analogy with a bore propagating upriver is 
very appealing, but further analysis is needed to see if the 


conditions necessary for a hydraulic jump can be verified. 


B. TEMPORAL VARIABILITY 

While the strength and occurrence of the internal bores 
show considerable temporal variability over the 34-d record 
considered here, there is no clear correspondence to the 
spring/neap cycle of the barotropic tide (Figure 19). The 
most regular bore passages occur during the neap phase at 
the end of the record, while an earlier neap phase does not 
show a clear bore signal. Clear bore signals are also seen 
during the first spring phase of the 34-d record while the 
second spring phase, with slightly smaller sea level ranges, 
has no clear bore signal. During the 01 to O09 September 
time period, when the bores are consistently strong and 
regular, they are phase locked (to within 1 h) to the sea 
level (Figure 18). In Figure 32, the 1026 kg/m isopycnal 
is taken to represent the background density field of the 


internal tide and the 1026.5 kg/m isopycnal is taken to 


Sa 


represent the isopycnal displacement of the internal bore. 
In these two cases, the background signal is phase locked to 
the sea level but the bore signal can be clearly seen out of 
phase with the background on 21 August (weak bottom 
intensification), and in phase on 09 August (strong bore). 
From this we can infer that strong bores will occur 
approximately 8.6 h after high tide when conditions are 


conducive to bore formation. 


Oe RECOMMENDATIONS AND IMPLICATIONS 

This thesis presents an initial description of some of 
the most notable features in the moored data from the 
Monterey Canyon internal wave experiment. There are still a 
number of unresolved issues however, including whether it 
can be determined to a high degree of certainty that what 
has been described here is in fact an internal bore. 
Secondly, while some intriguing connections between the 
background stratification and the strength and regularity of 
the bores have been noted, the dynamical reasons for this 
have yet to be worked out. Furthermore, the cause of the 
finger-like current pulses is unresolved. Also, the high 
frequency pulses that are moving down-canyon on 09 August 
have not been addressed in any detail. It is possible that 
they are internal waves formed by some topographic feature 
farther up-canyon, but the question remains as to why they 


are not seen every time there is down-canyon flow. Both of 


az 


these phenomena could be better explained by taking new 
measurements that have higher temporal and spatial 
resolution. An experiment that included continuous density 
measurements Ehueughouwe the water column = at several 
locations along the canyon axis would add insight into the 
development, propagation and evolution of the bores. 
Simultaneous sampling over a large section of the canyon 
would require numerous moorings and instruments, probably a 
cost-prohibitive exercise, but successive deployments in 
various areas farther down-canyon from the experiment site 
discussed in this thesis could accomplish many of the same 
goals. 

Even without a clear understanding of the formation and 
dynamics of the Monterey Canyon bores, the effects they have 
on measurements and operations in and around the canyon can 
be planned for. Several oceanographic institutions are 
located near Monterey Bay that conduct research in the 
canyon. One in particular, MBARI, uses ROV's to explore the 
canyon and is beginning to experiment with AUV's as well. 
Understanding and being able to predict the presence of 
strong internal bores can benefit these operations by 
reducing the risk of possible equipment loss. As recently 
as September 1999, an AUV was lost in the canyon possibly 
due to interference caused by an internal bore overtaking 
the vehicle. The effects of the bores can cause 


difficulties in operating in canyons beyond the threat of 


aS 


undesired motion. The bores are capable of resuspending a 
Significant amount of sediment as they propagate up the 
canyon (Bogucki et al. less SAF ase Significantly reducing 
alee a. tea Canyons also have a rectifying effect on the 
tidal currents in them (Lafuente et al. 1999) that can have 
the effect of bringing deep waters onto the shelf. Mixing 
associated with strong internal wave motions may also 
introduce high nutrient concentrations into the euphotic 


zone, thus enhancing productivity. 
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VI. CONCLUSIONS 

Thirty-four day records of near-bottom temperature and 
horizontal velocity throughout the lower third of the water 
column from within Monterey Canyon were examined. The 
internal tide was found to be highly non-linear with kinetic 
energy dispersed among numerous overtides near the bottom. 
Higher in the water column, energy iS more concentrated into 
the primary semi-diurnal constituent, M2. Strem@ilay boeeton-— 
intensified up-canyon currents were observed. The bottom 
currents and temperature varied in strength over the 34-d 
record, taking on the characteristics of an internal bore 
when they were at their strongest. Large up-canyon 
acceleration was accompanied by a rapid temperature drop. 
Figure 31 is a conceptualization of the internal bore that 
sometimes occurs at ~12-h intervals. 

When strong bores are present, they are nearly phase 
locked to the sea level variations and arrive at the 
measurement site about 8.6 h after the high tide at 
Monterey. They have maximum currents of ~55 cm/s with 
speeds exceeding 30 cm/s over a duration of approximately 
two hours. They are evident in the velocity record up to at 
least 35 m above the bottom, and may be accompanied by high 


frequency pulses extending even higher. Temperature drops 
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maximum acceleration (~0.24 cm/s’) signaling the bore 
arrival. 

The variation in the strength of the bottom-intensified 
currents is not caused by direct forcing from the tidal 
range like tidal bores on rivers (Pugh 1987), as the 
strongest bores do not occur during the spring phase of the 
spring/neap cycle. However, neither can it be said that the 
strongest bores will always be found during the neap phase 
of the cycle, even though that is the case during one 
POmeTon Of the records. Speculation on the reason for the 
strong variability observed in the strength of the bottom 
currents and the presence or absence of a bore centers on 
the stratification and the background shear in the mid-water 
eo /amn - The stratification on 09 and 21 August, days on 
which bores were present and absent, respectively, shows 
large differences. The mid- and deep water isopycnals 
oscillated in phase on the 9" and out of phase on the 21*, 
much as was observed by Petruncio et al. (1998) on two days 
in April and October 1994. Changes in the background 
stratification, perhaps caused by mesoscale events, could 
cause conditions more or less favorable to bore formation by 
changing the internal tide generation site and/or 
interactions with the background flow that would affect the 
strength of the internal wave and degree of wave steepening, 


a precursor to bore formation. 
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Figure 1. Visualization of internal wave propagation; phase 
and group velocities are along vectors k and C oF 
respectively (Munk 1981). 
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Figure 2. Illustration of internal wave reflection geometry 


(Slinn and Riley 1996), 9 here refers to the angle of energy 
propagation while in the text it refers to the angle of 


phase propagation. 
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Figure 3. Diagram of a 2-layer hydraulic jump (Holloway 
1987). This model can describe either a stationary ora 
propagating hydraulic jump. 
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Monterey Submarine Canyon 
Internal Wave Experiment 


August 1997 


= CID stations 


e ADCP moorings 


Depth contours in m 


BRS 


20 Kilometers 


l K Rosenjeld 





Figure 4. Experiment map. CTD station between the two 
moorings is referred to as "Shallow", other CTD station as 
iT] Deep Ww 
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Figure 5. High-resolution bathymetry data from MBARI 
showing locations of NB and BB relative to the canyon axis. 
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Figure 6. Concept for calculating principal axes for 
current data. +A' is considered up-canyon and -A" is 
considered down-canyon. 
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Figure 7. Current data from BB ADCP bin #1 (blue) with the 
derived principal axes from method one (green) and method 
two (red) overlain. 
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Figure 8. Comparison of BB ADCP bin #1 speed data (*) with 
method one speed (square) and method two speed (diamond). 
Where negative values indicate westward or down-canyon flow. 
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Figure 10. Same as Figure 9 except for BB ADCP. 
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Figure 11. Kinetic energy spectrum for each BB ADCP bin, 
with successive bins offset by one decade. The position of 
several tidal constituents and the Brunt-Vaisala frequency 
12 m above the bottom are indicated on the x-axis. 
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Figure 12. Kinetic energy spectrum from the near-bottom (12 
mab) and middle water column (124 mab) showing 
intensification compared to the GM79 open ocean internal 
wave kinetic energy spectrum. 
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Figure 13. 
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BB and NB temperatures with Monterey sea level 
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series. 
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Figure 14. Alpha Omega thermistor records, BB (black) and 
NB (gray) with Monterey sea level overlain for 09 August 
1997. 
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Figure 16. BB ADCP pitch sensor (X), along-canyon velocity 
(squares), and Alpha Omega thermistor (*) on 09 August 1997 
showing bore passage. 
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Figure 17. Differences between successive high tides and 
internal bores (01 to 09 September 1997). 
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Figure 18. Time delay of internal bore passage relative to 
previous high tide. 
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Figure 19. Sea Level ranges (black) and BB ADCP pitch 
sensor (gray). 
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Figure 21. Speed contours (color, cm/s) from BB mooring and 
VM ADCP, O09 August 1997, showing the three events described 
in ch. IV with isopycnal contours (black lines, .025 kg/m 
spacing, min 1026 / max 1026.675 kg/m) from CTD casts 
conducted at Shallow CTD station and temperature from the 
Alpha Omega thermistor (heavy black). 


qa. 





E/W velocity (cm/s) 





00:00 06:00 12:00 18:00 00:00 


N/S velocity (cvs) 





00:00 06:00 12:00 18:00 00:00 


Figure 22. Comparison of VM-ADCP 1-min (squares) at 215 m, 
VM~-ADCP 3-min (diamonds) at 215 m, and BB-ADCP (*) at 217 m. 
Upper(Lower) panel is E/W(N/S) component of velocity. 
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Figure 23. BB ADCP horizontal current speed (vector 
magnitude) and direction (true north up) for 09 August 1997. 
Maximum speed plotted is 69 cm/s. 
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Figure 24. Speed contours (color, cm/s) from BB ADCP, 09 
August 1997, showing the 0700 event with isopycnal contours 
(black, 0.025 kg/m spacing, min 1026.025 kg/m / max 
1026.675 kg/m) from CTD casts conducted at Shallow CTD 
station and temperature from the Alpha Omega thermistor 
(heavy black). 
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Figure 25. Same as Figure 24 showing the 1845 event on 09 
August 1997 in detail. Minimum density plotted is 1026 
kg/m; maximum density plotted is 1026.475 kg/m. 
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Figure 26. 


min period) down-canyon current pulses on 09 August 1997 in 


detail. 
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Figure 27. 
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Figure 28. Same as Figure 24 showing 21 August 1997. 
Maximum density is 1026.675 kg/m; minimum density is 
1026.200 kg/m’. 
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Figure 30. BB ADCP horizontal current speed (vector 
magnitude) and direction (true north up) for 21 August 1997. 
Maximum current plotted is 44 cm/s. 
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Figure 31. Conceptual image of the internal bore. T1 
through T4 represent isotherms (0.25 °C apart) which 
terminate ina front in the up canyon direction. Lower 
frame is blow-up of Ti in the region of the bore showing 
turbulence-induced structure. 
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Figure 32. Sea level and isopycnal depths for 26 and 26.5 
Sigma t surfaces on 09 (black) and 21 August (gray). 
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